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IS QUASAR VARIABILITY A DAMPED RANDOM WALK? 
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ABSTRACT 

The damped random walk (DRW) model is increasingly used to model the variability in quasar light 
curves, but it is still uncertain whether the DRW model provides an adequate description of quasar 
variability across all time scales. Using a sample of OGLE quasar light curves, we consider four 
modifications to the DRW model by introducing additional parameters into the covariance function 
to search for deviations from the DRW model on both short and long time scales. We find excellent 
agreement with the DRW model on time scales that are well sampled by the data (from a month to 
a few years), possibly with some intrinsic scatter in the additional parameters. On very short time 
scales (below a few months), we see some evidence of the existence of a cutof f time scale below which 
the correlation is stronger than the DRW model, echoing the recent finding of Mushot zkv et al.l (|2011[ ) 
using quasar light curves from Kepler. On very long time scales (> a few years), the light curves do 
not constrain models well, but are generally consistent with the DRW model. 

Subject headings: galaxies: active — galaxies: statistics — methods: data analysis — methods: 
numerical — methods: statistical 



1. INTRODUCTION 

The optical variability of quasars has long been a pro- 
posed diagnostic of the central engines of quasars (e.g., 
iKawaguchi et al.l I1998D . While the physical origins 
of variability rema ins an open question (see, e.g., 
iDexter fc Agol|[20Tlh , recently it has been proposed that 
the optical variability of quasars is mathem atically well 
descr ibed by a damped random walk (DRW, iKellv et all 
|2009f ). which characterizes quasar light curves as a 
stochastic process with an exponential covariance func- 
tion S(At) = a 2 exp(—\At/r\), defin ed by an amplitud e 



a and a charact e ristic time scale r. IKellv et al.l ( 2009h . 
iKozlowski et al.l (|2010t l and iMacLeod et al.l (|201()h have 



shown that the DRW model provides a viable expla- 
nation for the variability of individual quasars found 
in a heterogeneous sample, OGLE (Optical Gr avita- 
tional Lcnsing Exp eriment. lUdalski et all |l997) and 
SDSS Stripe 82 (S82. lSesar et al.ll2007f ). respectively, and 
that the parameters r and a are correlated with the phys- 
ical properties of the quas ar (rest wavelength, lum inos- 
ity, and black hole mass). IMacLeod et al.l (|2011aD also 
demonstrated that the ensemble structure functions of 
quasars can be recovered by modeling each individual 
quasar with a DRW based on its luminosity, rest wave- 
length and estimated black hole mass. 

The DRW model is also a powerful tool. For exam- 
ple, it provides a well defined approach to variability- 
selecting quasars (IKozlowski et alj|2010t IMacLeod et al.l 
I2011bt iButler fc Bloom! 120 111) . It also provides a well- 
defined approach to interpolating and modeli ng quasar 
light c urves in reverberation mapping studies. IZu et all 
(|2011h show how it provides a general approach to mea- 

Department of Astronomy, The Ohio State Univer- 
sity, 140 West 18th Avenue, Columbus, OH 43210, USA; 
yingzu@astronomy.ohio-state.edu 

2 The Center for Cosmology and Astroparticle Physics, The 
Ohio State University, 191 West Woodruff Avenue, Columbus, 
OH 43210, USA 

3 Warsaw University Observatory, Al. Ujazdowskic 4. 00-478 
Warszawa, Poland 



suring lags in one or more emission lines or velocity bins 
of a si ngle line. The model has b een used bv iGrier et al.1 
(|2012f ). and a variant was used bv lPancoast et al.1 (|2011l ). 

However, it is still uncertain whether the DRW model, 
with only two parameters, can fully characterize quasar 
variability across all ti mescales. Early studies (e.g. 
ICollier fc Peterson! 1200 1|) examined the power spectrum 
density (PSD) of a small number of individual quasars, 
finding evidence for the f~ 2 power law behavior of 
the DRW model on short time scales, and weaker evi- 
dence for a flattening on long time sca l es. The first se- 
ries of DRW studies (IKellv et al.l [2009t IKozlowski et al.l 
120101 : IMacLeod et al.ll2010l ) essentially followed these re- 
sults and reconfirmed the average PSD structure, but 
did not explore any possible deviations. We know 



from lMushotzkv et al.l (|2011l ) that the DRW model prob- 
ably fails on very short time scales (< week), hav- 
ing too much power on short time scales compared to 
that measured by Kepler. Mathematically speaking, the 
DRW model (a.k.a. the Ornstein-Uhlenbeck process) 
is unique — it is the only statio nary Gaussi a n pro cess 
that is "memoryless" (Markov). IKellv et all (|2009f ) in- 
troduced DRW model as a first order continuous au- 
toregressive (CAR(l)) process, and its discrete counter- 
part (AR(1)) is simply an iteration Xi+i = aARfi)^; +£j, 
where e is a Gaussian deviate. Physically, if quasar vari- 
ability is driven by multiple mechanisms, or multiple 
seeds of the same mechanism with different relaxation 
time scales or amplitudes, there should be deviations 
from the DRW model. 

In this paper, we explore several alternative covariance 
functions with a third parameter that allows deviations 
from the DRW model and explore whether they provide a 
better representation of quasar variability. This needs to 
be done for individual quasars because we are now certain 
that the ensemble statistics (e.g., structure functions) of 
quasars are weighted averages o f individual quasars wit h 
different individual statistics (IMacLeod et al.l l2011bD . 
We employ the OGLE light curves of quasars behind 
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the Small and the Large Magellanic Cloud (SMC and 
LMC). Thanks to the high-cadence and long-baseline of 
OGLE, the light curves generally have ~ 570 photomet- 
ric epochs over ~ 7 yrs. making them extremely well- 
suited for studies of quasar variability. We describe this 
data set and our sample selection in Section [2j In Sec- 
tion [3] we introduce four new covariance functions that 
we use to test for deviations from the DRW model and 
describe the model fitting procedures. We present our 
main results in Section |U We conclude by discussing the 
physical implications in Section [5] 

2. OGLE QUASAR LIGHT CURVES 

We use the light curves of q uasars behind the LMC an d 
SMC monitored by OGLE (lUdalski et al.l fl997l [2008T) . 
Most of these were identified by iKozlowski et al.1 (|201lL 
l2012h in part fro m candidates variabili ty-selected using 
the DRW model ijKozlowski et al.ll2010f ). The photomet- 
ric uncertainties of the light curves we re estimated b y 
using Difference Image Analysis (DIA, lWozniak|[2000l ). 
The initial sample consists of 223 I-band quasar light 
curves. There are typically ~ 570 epochs taken over ~ 7 
years on a 2 day cadence with six month seasonal gaps 
when OGLE instead focuses on the Galactic bulge. For 
comparison, the S82 quasars have only ~ 50 epochs nom- 
inally over ~ 10 years, but there are really a few epochs 
in 1998, a gap until 2000 with 1-3 epochs per year cover- 
age to 2005, and then 10-30 epochs per year from 2005 to 
2008. The OGLE quasars should allow us to investigate 
quasar variability on both small and long timescales in 
more detail than the SDSS S82 light curves. 

Since we are examining the problem of additional pa- 
rameters, we need higher quality light curves than if we 
were only trying to estimate DRW parameters. Start- 
ing from the 223 light curves, we first selected sources 
with high S /N for their variability, keeping sources where 
the rms magnitude variation cr rms was significantly larger 
than the mean photometric error a m (cr lnls /a m > 2.0). 
We also required a m < 0.1 mag. This left us with 87 light 
curves in the sample. Next we required that the light 
curves yield a well-defined DRW time-scale r by requiring 
that both ln£o/£ max < —0.5 and ln£ 00 /£ max < —0.5, 
where Co, and £ max are the likelihoods for r = 0, 
r = oo and the best-fit r, respectively, where the best- 
fit t are in the range of 17 < r < 2700 days. This is to 
ensure that the light curves have some statistically de- 
tectable characteristic timescale. The final sample has 55 
light curves of quasars with redshifts between 0.15 and 
2.50 and I-band apparent magnitudes between 16.7 and 
19.4 mag. 

About 45% of the objects are at redshifts where the 
I-band continuum variabilities are contaminated by the 
lagged variabilities of broad emission lines. For a broad 
band filter like I-band and assum ing a typical quasar 
spectrum (|Vanden Berk et al"1l2001[ ). the line contributes 
only < 6% of the flux in the filter (except for two objects 
at z ~ 2.45 where the Ha line contributes upto 14% 
of the I-band flux). Since we are probing the shape of 
the covariance functions instead of the variability ampli- 
tudes, we do not expect this small level of contamination 
to affect our results. We tested for any effects by divid- 
ing our sample into line-contaminated and line-free sub- 
samples, and found negligible differences between the two 
sub-samples in the likelihood analyses. Since the results 
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Figure 1. Comparison of the powered-exponcntial (panel a), 
Matcrn (panel 6), Pareto-exponential (panel c), and Kepler- 
exponential (panel d) covariance functions (thin solid curves) with 
the DRW model (thick dashed curves). The third parameter of 
each model is given by the gray-scale colorbar in each panel. The 
models arc normalized to unity at zero lag. 

seem little affected by broad line contamination, we will 
not discuss this issue further. 

3. COVARIANCE FUNCTIONS AND ANALYSIS METHOD 

We model each light curve as a Gaussian process, de- 
fined by a covariance function S(At) = cov(m(t),m(t + 
At)) where At is the time separating two observations 
and we assume a stationary process with a well-defined 
mean. The covariance function is related to the structure 
function (SF) by 



SF(At) = y/2a 2 - 2S(At) 
and to the PSD by a Fourier Transform 



P(f) = J S(At)e~ 2 * Atfl dAt. 



(1) 



(2) 



In particular, the DRW model has a covariance func- 
tion of 

S mw (At) = a 2 exp ( -] (3) 



and a PSD of 



P(.f) = 



4a t 



1 



(27TT/) 5 



(4) 



We do not repeat the SFs of the individual models be- 
cause they are trivially related to the covariance func- 
tions (Equation [T|). The DRW model can be equiva- 
lcntly parameterized by r and the asymptotic amplitude 
of the SF (SFoo = V2cr, IMacLeod et aLll2?)Tnh . or r and 
the slope of the SF on short time scales (a = \/2a 2 /r, 
IKellv et al.l[2009l ). 

We consider the four modified covariance func- 
tions shown in Figure [1] A natural variant of the 
DRW covariance function is the so-c alled "powered- 
exponential" (PE) function (|Neallll997t) 



S PE (At) = a 2 exp 



dt 



for < 7 < 2, (5) 
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Figure 2. Combined (left panels: a and c) and individual (right 
panels: b and d) marginal log-likelihoods for 7 in the powered- 
cxponential model for the data (top panels: a and b) and the mock 
data (bottom panels: c and d). For the DRW model 7 = 1.0, as 
indicated by the vertical dashed line. 

where the special case 7 = 1 corresponds to the DRW 
model. The PSD for the PE covariance can be ex- 
pressed in a closed form only for 7 = 1 (DRW) and 
2 (Gaussian), but for small time scales (s — > 00), P(s) oc 
s-( 7+1 ) dLim fc Muniandvll2003D . The PE model modi- 
fies the DRW model on both long and sm all scales (panel 
a of F igure Q]) and is the model used in iPancoast et al.1 

(EqU. , 

Th e "Matcrn" (MA) covariance function (jMaternl 
HflMl is defined by 



'MA 



(At) 



Figure 3. The effects of incorrect error estimates for the powered— 
exponential model. The top (bottom) panels assume that the true 
photometric errors are over (under) estimated by 20%. As in Fig- 
ure [2] the left (right) panels show the joint (individual) estimate 
of the new parameter 7. 



for a short time scale cutoff. To search for the cutoff 
within our sample light curves, we define the Kepler- 
exponcntial (KE) covariance function by modifying the 
DRW model on very short time scales to be 



5 KE (Ai) = a 2 



1 - 
exp 



(1 — exp 
(-1^1) 



1^1)) 



a/ 



3/2 



At < T cut 
At > r cut , 

(8) 



with v 



(6) 

where T is the gamma function, K v is a modi- 
fied Bessel function of the second kind, and the 
new parameter, v primarily adjusts the correlation 
function on short time scales. The covariance func- 
tion is exponential on long time scales SMA(At) = 

u 2 ^Y(v)- 1 ( s /vJ%At\lry- 1 ^%^-^^ t \l T as 
At — » 00. The MA covaria nce function has a 
PSD (jAbramowitz fc SteguriflQfiSh of 



where r cu t is the cutoff time scale below which the covari- 
ance function is stronger than the DRW model, so that 
> ihc PSD has a similarly steep slope as found for the Ke- 
ller light curves. The model returns to the DRW model 
for r cut = (panel d of Figure [T|) , although we should 
not be able to recognize values of r cu t on scales smaller 
than a few times the typical cadence (~ 10 days). 

In order to test for deviations from the DRW on 
long timescales, we combined the DR W model on sh ort 
timescales with the "Pareto" function ()Paretol[l967t ) on 
long timescales (the "Pareto-exponential" , PA), 



S PA (At) = a 2 



P(f) 



2tt tT(v + 1/2) 



2i> 



exp( 

1 I At |- Q 



A/ 



At < T 
At > T 



(9) 



(7) 



The special case v = 0.5 gives the DRW covariance func- 
ti on (panel b of Fig u re [It. 

iMushotzkv et al.l (|2011[ ) found that the quasar light 
curves measured by Kepler do deviate from the DRW 
model on very short time scales (~ days), with steeper 
PSD power-law slopes of —2.6 to —3.3, rather than the 
-2 of the DRW model (Equation HJ). Thus, there must 
be a T cu t below which the fluctuation amplitudes cut 
off more sharply than in the DRW model. The OGLE 
light curves, with a typical cadence of ~ 2 days, over- 
lap the time scales of the Kepler light curves (~ 6 hr 
to ~ 1 month), so there is some prospect of searching 



where the power-law index a > 1 controls the extent of 
the covariance tails (panel c of Figure [T]). In this case, 
no value of a exactly matches the DRW and placing the 
break at time scale r is reasonable but arbitrary. 

Our goal is simply to estimate a value for the new 
parameters introduced in each of these models using a 
Maximum Likelihood (ML) approach. We briefly recap 
the compo nents of the like li hood function he r e and refer 
readers to IZu et all (120111) . iKozlowski et al.l (|2010D and 
IRvbicki fc Press! (|1992fl for details. We model each light 
curve as 

m(t) = s(t) + n + Lq, (10) 
where s(t) is the underlying variability signal with covari- 
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Figure 4. As in Figurc[2]but for the Matern model. Here v = 0.5 
corresponds to the DRW model. 

ancc matrix 50, n is the measurement uncertainty with 
covariance matrix iV, L is a vector with all elements equal 
to one, and q is the light curve mear0. After optimizing 
the value of q, the likelihood of the model parameters is 



£(t,ct, [ 7 ,^,r cut ,a]) = iq-^I^C-^r^expf - 



m T CT 



where C is the overall data covariance S + N and 
C- 1 - C- 1 L(L T C- 1 L)- 1 L T C- 1 . 



ci 1 



(11) 



(12) 



We first fit each quasar light curve with the DRW 
model to obtain Maximum Likelihood (ML) estimates 
for a and r, which are then used to generate a mock light 
curve that is fully consistent with the DRW model while 
having exactly the same sampling and noise properties 
as the true light curve. We use these as a comparison 
sample when evaluating the evidence for any additional 
parameters. Next, we fit both the data and the mock 
light curves using each of the four covariance functions. 
Since we only care about the differences in the third pa- 
rameter, we use a Maximum Likelihood analysis based 
on £([7, [7, 1/, r cut , a]) where 

Tmax and Oraax maximize the likelihood at fixed 3rd pa- 
rameter for each light curve. 

4. RESULTS 

The results of our model depend on how well the 
DIA photometric errors of the light curves in our sam- 
ple are determined. Assuming the DRW is an adequate 
model for each individual light curv^B, we can estimate 
the fractional errors in the photometric error estimates 
from the % 2 /dof distribution of the DRW fits to the 55 

4 The entries of Sij are simply the values of the covariance func- 
tion Sij = S(Atij), so we have used the same symbol for both 

5 The simultaneous fit of q is important because any con- 
stant level of contamination (e.g., host galaxy light) is removed 
by marginalizing over q. 

6 The DRW model has one fewer degree of freedom than other 
models, so the inferred error statistics should be more conservative 
than using other models 



Figure 5. Distribution of quasars as a function of the best— fit 
Tcut for the data (left panel) and mock (right panel) light curves 
using the Kepler-exponential model. The vertical dashed line in 
each panel at 10 days indicates the scale below which no variabil- 
ity signal can be resolved for the lack of sampling the Kepler- 
exponential covariance function below r cu t- 

light curves. For all these objects, the overall good- 
ness of fits to the light curves with the DRW model 
are reasonable (0.88 < x 2 /dof < 1.53), with an average 
{x 2 /dof) ~ 1.09 and dispersion a x 2 / do f = 0.13, while we 

would expect (y/ 2 /dof) = 1 and a x 2 j do f = ^2/ dof ~ 
0.05. These differences could be explained by fractional 
shifts a m = er^" e (l + e) between the true uncertainties 
a^ ue and the nominal error bar a m where the mean frac- 
tional error is (e) = —0.04 with a dispersion of a e = 0.06. 
This will help us understand the scatter of the additional 
parameter between quasars. 
L We simply discuss each of the four models in turn, 
^Starting with the PE model results shown in Figure [2j 
The panels on the left show that the combined results 
agree with the DRW model remarkably well. The joint 
marginal likelihood function for 7 strongly peaks at 1.0 
with ( 7data ) = 0.98 ± 0.04 using Aln(£/£ max ) = -2.0 
for the error estimate (which would correspond to 2a 
for Gaussian uncertainties). This is generally consistent 
with the result for the mock light curves of (7mock) = 
1.04 ± 0.03. The individual likelihoods, shown in the 
right panels of Figure El seem to show a larger scat- 
ter in the data (cr 7l data — 0.37) than in the mock 
data (cr 7 , m ock = 0.19). While essentially all the light 
curves are consistent with the DRW at better than 
3cr ((A ln(£/£ max ) = —4.2)), there is one extreme out- 
lier (seen as the spike at 7 = 0.14). While there is noth- 
ing obviously odd in this object's light curve, it is very 
close to our noise selection limit and a source will be 
biased towards low 7 (close to white noise) if the true 
measurement errors are larger than estimated (see be- 
low). This object is also an outlier in the other models. 
Nonetheless, excluding or including it has a negligible 
effect on the combined likelihood functions. 

Before going to other models, we want to interpret the 
larger spread we observe in the best-fit 7 for the individ- 
ual data light curves (panel b of Figure [2]) . The differ- 
ence of the dispersions between the individual best-fits 
in the data and mock light curves is 07 = (a 2 data — 

07, mock) 1 / 2 = 0-31 (0.29 if we drop the one extreme out- 
lier). This can be interpreted either as evidence for in- 
trinsic scatter or as evidence for problems in the pho- 
tometric uncertainties a m . As an experiment, we re-fit 
the PE model to the mock light curves after resetting 
a m to be 20% larger (e = +0.2) or smaller (e = —0.2) 
than the actual uncertainties used to generate the mock 
light curves. The results of the experiment are shown in 
Figure [3l When the noise is over-estimated, true vari- 
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Figure 6. As in Figure P21 but for the Pareto-exponcntial model. 
No value of /3 corresponds to the DRW model. 

ability power on short time scales is instead interpreted 
as noise, so the models shift to higher 7 for stronger cor- 
relations on short time scales. When the noise is under- 
estimated, noise on short time scales is interpreted as 
signal, so the models shift to lower 7 for weaker correla- 
tions. As expected, in Figure El over-estimating the er- 
rors biases the estimates of 7 to be high, with (7^ ock ) = 
1.30 ± 0.04 (panel a of Figure [3]), and vice versa when 
under-estimating errors, (7~ ock ) = 0.62 ± 0.01 (panel c 
of Figured]). Given that a 20% error in a m causes a 0.3 
bias in the estimate of 7, so that 7 ~ 1 + 1.5e, then the 
excess variance in the estimates of 7 for the data could 
be explained by making o~ e = 0.20, which is significantly 
larger than the a e = 0.06 suggested by the x /do/ dis- 
tribution. Also note that the estimate of (e) from the 
\ 2 /dof distribution would tend to shift the estimates of 
7 for a DRW model to 7 = 1 + 1.5 (e) ~ 0.93, so sys- 
tematic uncertainties in the photometric errors produce 
uncertainties comparable to the statistical uncertainties 
in (7). If we use o~ e = 0.06 from the \ 2 /dof distribution 
as an estimate of the scatter in the fractional uncertain- 
ties in a m , then we appear to be left with a'" 1 = 0.28 of 
intrinsic scatter between quasars. 

Similarly, Figure 2] shows the results for the MA model. 
The likelihood function for the data is again strongly 
peaked at the value that corresponds to the DRW model, 
v = 0.5, with (t'data) = 0.492 ± 0.016, as compared to 
(iVock) = 0.525 ± 0.017. The distribution of the indi- 
vidual results, is again broader for the real data than for 
the mock data, suggesting an excess scatter of a v = 0.08 
in v if there are no other systematic uncertainties. In 
this case, the sensitivity to the uncertainties in the pho- 
tometric errors is v = 0.5 + l.le, suggesting an intrinsic 
scatter of 4" = 0.05. 

On even shorter time scales, the histograms in Fig- 
ure[5]show the distributions of the best-fit r cut of the KE 
model for the data (left panel) and the mock (right panel) 
light curves. For the data, there is a clearly bimodal dis- 
tribution of r cut demarcated by a gap around 10 days 
that corresponds to 5 times the typical cadence (2 days). 
The first peak below 10 days comprises objects that are 



either consistent with the DRW model (r cut = 0) or have 
a r cut that is unresolved by the OGLE sampling (^ 4 
discrete values of At below r cut could be computed for 
Ske in Equation[8|). The second peak at 30 — 100 days is 
well beyond the cadence scale, indicating another popu- 
lation of objects (23 out of 55) that favor a cutoff time 
scale beyond a month. For the mock sample, there is 
only a single peak below 10 days and a steady decrease 
towards longer T cut , with only 11 out of 55 objects having 
r cu t larger than a month. Unfortunately, the estimates of 
r cu t are very susceptible to errors in the estimated photo- 
metric errors, so the second peak in the r cut distribution 
of the data light curves could also be an artifact of the 
scatter in e. Thus a cutoff time scale r cut ~ 1-3 months 
may be marginally detected in half of the quasars, but 
the cadence and photometric uncertainties of the data 
are not optimal for probing this regime. 

Finally, the light curves do not constrain changes in 
the structure of covariance function on longer timescales. 
The individual likelihood functions for a in the Pareto- 
exponcntial model are mostly flat, as shown in Figure [6j 
The combined constraint from the data favors a slightly 
higher value of (adata) = 3.38 ± 0.45 than the mock light 
curves, (a moc k) = 2.83 ± 0.42, but the two peaks are 
mutually consistent. The uncertainties are, however, too 
large to justify additional models for the long time scale 
behaviors. 

5. DISCUSSION 

We considered four families of covariance functions to 
test whether there is any evidence that quasar variability 
differs from a DRW model on time scales of weeks to 
years. Using a set of well-sampled OGLE quasar light 
curves, we find no indication of a third parameter in three 
cases (PE, MA, and PA), but evidence for the existence 
of a cutoff time scale in the KE case. More specifically, 

• On the time scales best sampled by the light 
curves (from months to years), the typical light 
curve is clearly well-described by the DRW model, 
potentially with some scatter in the true structure 
of the power spectrum if the photometric error es- 
timates are correct to better than \o- m /a t ^ l ue — 1| < 
10%. 

• On very short time scales, there are hints of a char- 
acteristic cutoff time scale r cut (1-3 months) for 
~ 50% of the quasars, below which the correla- 
tions become stronger than predicted by the DRW 
model. The OGLE data a re clearly not comp etitive 
with the Kepler study by iMushotzkv et al.l (|2011l ) 
in this regime. 

• On very long time scales, the light curves are con- 
sistent with the DRW model, but the precision of 
the test is limited by the length (7 years) of the 
light curves. 

The apparent scatter about the DRW model may well 
be evidence for multiple stochastic processes rather than 
a single DRW. For example, the powered-cxponential co- 
variance function can be viewed as the mixture of a con- 
tinuous sum of independent DRW models with charac- 
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where s = tj/t, (7 2 (st) is the amplitude of process i, and 
P(s, 7) is the probability density function for the pro- 
cess on time scales r. The combination of P(s, r y)a 2 (sT) 
is then the inverse Laplace transform of the powered ex- 
ponential (|Johnstodl2006h . For 7 = 1 (DRW), P(s, 7 ) 
is a Dirac S function at s = 1, so there is only one pro- 
cess operating at one timescale r. Adding several weaker 
processes with different time scales would then appear as 
a scatter about 7=1, but would be difficult to distin- 
guish from a sin gle DRW model. T his is similar in spirit 
with the work of iKellv et al.1 (|2011l ) , where they modeled 
quasar optical variability as a linear mixture of different 
DRW processes and found it generally provides no better 
fit t han a single DRW mo del. In the toy model adopted 
by iDexter fc Agoll (|2011l ). they matched the observed 
DRW variability by exciting many independent fluctu- 
ating zones with exactly the same r = 200 days in the 
accretion disc, and then argued that the resulting change 
in the effective temperature profile of the disk would ex- 
plain the observed discrepancies between thin disc sizes 
inferred from gravitational mic rolensing of lensed qu asars 
and their optical luminosities (| Morgan et al.ll2010h . 

On short time scales, it would be useful to have oc- 
casional campaigns in which OGLE sampled shorter 
time scales with longer integration times to fill in these 
time baselines with higher precision data. On long 
time scales, the OGLE-IV project is already extend- 
ing the light curves, so it will become increasingly pos- 
sible to explore the long time scale behavior of indi- 
vi dual quasars. T h e con tinuing spectroscopic surveys 
by iKozlowski et all (|2012[ ) will also greatly increase the 
numb er of confirmed qua sars wit h OGLE light cu rves. 
DES (ICollaborationl[2005h. LSST (llvezic et al.l[200l and 
Pan-STARRS (jKaiser et all 120021 ) will carry out similar 
extensions for the SDSS Stripe 82 quasars but generally 
at lower cadence (albeit for larger numbers of objects). 
Given a larger sample it wll also be possible to search for 
correlations of the apparent deviations from the DRW 
model with other physical characteristics of the quasars 
(luminosity, wavelength, etc.). 
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